PASJ: Publ. Astron. Soc. Japan , 1-??, - 
(gj 2011. Astronomical Society of Japan. 



Suzaku Metal Abundance Patterns in the Outflow Region of M82 and 

the Importance of Charge Exchange 



O 

(N 

< 
(N 

o 

u 

43 

9 L1, 
6 



> 
oo 

in 
o 



Saori Konami, 1 ' 2 Kyoko Matsushita, 1 Takeshi Go Tsuru, 3 Poshak Gandhi, 4 

and Toru Tamagawa 2 ' 1 
Department of Physics, Tokyo University of Science, 1-3 Kagurazaka, Shinjuku-ku, Tokyo 162-8601 
2 High Energy Astrophysics Laboratory, the Institute of Physical and Chemical Research, 2-1 Hirosawa, Wako, Saitama 

351-0198 
konamiQcrab.riken.jp 

3 Department of Physics, Kyoto University, Kitashirakawa-oiwake-cho, Sakyo-ku, Kyoto 606-8502 
4 Institute of Space and Astronautical Science (ISAS), Japan Aerospace Exploration Agency, 3-1-1 Yoshinodai, 



Chuo-ku, 



Sagamihara, Kanagawa 252-5210 

(Received ; accepted ) 
Abstract 

We performed spectral analysis of Suzaku data of the galactic disk and outflow regions of the starburst 
galaxy M82. Thermal modeling of the central disk regions requires at least three temperature components. 
The Ly/3 line fluxes of OVIII and NeX exceed those expected from a plasma in collisional ionization 
equilibrium. The ratios of Ly/3 /Lya lines for O VIII and Ne X are higher than those of collisional ionization 
equilibrium, which may be caused by the process of charge exchange. In the outflow wind region, the spectra 
are well reproduced with two-temperature thermal models, and we have derived the metal abundances of 0, 
Ne, Mg, and Fe in the outflow. The ratios of O/Fe, Ne/Fe, and Mg/Fe are about 2, 3, and 2, respectively, 
relative to the solar value determined by Lodders (2003). Since there is no evidence of charge exchange 
in outflow region, the metal abundances should be more reliable than those in the central region. This 
abundance pattern indicates that starburst activity enriches the outflow through SN II metal ejection into 
intergalactic space. 

Key words: galaxies: individual (M 82), galaxies: interstellar medium, galaxies: abundances, galaxies: 
starburst, X-rays: galaxies 



1. Introduction 

Metals in hot gas in the interstellar medium (ISM) 
and intergalactic medium can provide a cumulative fos- 
sil record on the history of star formation and evolution 
of galaxies. Starburst galaxies play an important role in 
supplying intergalactic metals via outflows heated by type 
II supernovae (SN II). These outflows are thought to con- 
tain metals synthesized by SN II and may have enriched 
the intra-cluster medium (ICM) and warm-hot intergalac- 
tic medium (WHIM). 

Recent X-ray observations of starburst galaxies have re- 
vealed extended halos of hot plasma around several galax- 
ies (e.g. Wang et al. 2005; Strickland 2007; Wang et 
al. 2009). However, measurements of metal abundances 
within these outflows are limited. For NGC 253, XMM- 
Newton detected the metals O, Ne, Mg, Si, and Fe in the 
central regions <3 kpc along the outflow axis (e.g. Bauer 
et al. 2007). For NGC 4631 the abundance patterns for 
O, Ne, Mg, Si and Fe were observed in the extended halo 
region, ^13 kpc with Suzaku (Yamasaki et al. 2009). The 
outflow wind and disk regions of NGC 4631 show pat- 
terns close to those expected from the SN II yields and 
solar abundance (Lodders 2003), respectively. 

M 82 is the most famous, prototypical starburst galaxy 
with an X-ray emitting outflow (e.g. Tsuru et al. 1997; 



Stevens et al. 2003; Origlia et al. 2004). Tsuru et al. (1997) 
observed M 82 with ASCA and derived a peculiar abun- 
dance pattern of O/Ne/Mg/Si/S and Fe in the hot X-ray 
emitting plasma around the starburst region. No combi- 
nation of SN la and SN II could reproduce the observed 
abundance pattern. Furthermore, Ranalli et al. (2008) 
reported spatial variations of the abundances with XMM- 
Newton. They found the abundances of O, Ne, Mg, Si, 
and Fe to increase from the center going out to the exter- 
nal outflow region, the absolute value of the O abundance 
they found was extremely low. In contrast, Tsuru et al. 
(2007) derived SN II like abundance pattern in the 'cap' 
region, which may be a large photoionized cloud in the 
halo of M82, at a distance of 11.6 kpc north from the 
center. 

Lallement (2004) has pointed out that line emission in 
the 'cap' region could arise from the process of charge- 
exchange (CX) , where highly ionized gas may be colliding 
with the cold cloud. Tsuru et al. (2007) reported that an 
emission line consistent with the C VI transition of n = 4 to 
1 at 0.459 keV has been marginally detected, although it 
is not statistically significant at the 99% confidence level; 
the presence of this line would suggest charge-exchange 
processes in the 'cap' region with Suzaku. In the XMM 
spectra in the central region, Ranalli et al. (2008) reported 
line-like features around ~ 0.7 keV and 1.2 keV, which 



2 



S. Konami et al. 



were interpreted as CX reactions of neutral Mg and Si, 
respectively. Most recently, in the XMM reflection grating 
spectra (RGS), Liu et al. (2011) reported detection of the 
CX emission of OVII, NelX, and MgXI triplets at the 
center region. An important piece of evidence of CX is 
that the forbidden lines are stronger than the resonance 
lines, especially for OVII ions. They concluded that a 
significant fraction of these lines originate as a result of 
CX emission rather than thermal radiation. 

CX occurs between highly ionized atoms and neutrals, 
and emits X-ray emission lines. When a highly ionized 
ion hits a neutral atom, an electron is captured in an ex- 
cited state in the ion. It then decays to the ground state, 
emitting X-ray photons. CX emission has been detected 
from regions where hot and cool gas co-exist. The first 
discovery of CX was from comet Hyakutake (Lisse et al. 
1996), following which X-ray emission was detected from 
many comets. More recently, the diffuse soft X-ray back- 
ground was found to contain CX emission from neutrals in 
the geocorana/heliosphere and solar wind ions (e.g. Cox 
1998; Cravens 2000). In studies of supernova remnants, 
CX emission can explain anomalously high intensity ra- 
tios (i.e. f3/a and 7/a lines) compared with those of col- 
lisional ionization equilibrium (CIE) plasma (Rasmussen 
et al. 2001, Katsuda et al. 2011). In the M 31 bulge, an 
intensity excess at the O VII triplet was observed owing to 
CX emission (Liu et al. 2010). 

Over the 3 keV energy range in the spectra of M 82, 
there are hard components whose nature is not yet fully 
understood. Some part of this emission comes from 
point sources (X-l, for example; Miyawaki et al. 2009). 
With Chandra and XMM-Ncwton, Strickland & Heckman 
(2007) detected diffuse hard X-ray emission and Fe-K line 
emission which cannot be attributed to CIE plasma. 

In this paper, we investigated the spatial distribution 
of abundance patterns in the outflow wind of M 82, be- 
tween the center and the 'cap' regions. We also derived 
ratios of K/3 to Ka lines of O and Ne and studied the 
spatial variation of the effect of CX emission. We expect 
both hot ions and cold atoms exist in the center and the 
'cap' regions. This paper is structured as follows. In sec- 
tion 2, we summarize observations of Suzaku and XMM- 
Newton. Sections 3 and 4 detail the data analysis and re- 
sults. Section 5 gives a discussion of these results. Finally, 
in section 6 we present our conclusions. 

We adopted 3.53 Mpc for the distance to M 82, which 
was determined by the luminosity of the tip of red giant 
branch (Karachentsev et al. 2004). Unless noted other- 
wise, we use the new solar abundances in Lodders (2003), 
and the quoted errors are for a 90% confidence interval 
for a single interesting parameter. 

2. Observation and Data Reduction 

2.1. Suzaku 

Suzaku observed M 82 on three occasions in 2005 
October during the Science Working Group phase to in- 
vestigate the outflowing hot gas from the galactic plane. 
The logs of three observations are listed in tablel. The 



XIS consists of three front-illuminated (FI: XIS0, XIS2 
and XIS3) CCD cameras and one back-illuminated (BI: 
XIS1) CCD camera (Koyama et al. 2007). Combined with 
XRT (Serlcmitsos et al. 2007), the field of view (FOV) of 
the XIS covers a ^18'xl8' region with a half-power di- 
ameter (HPD) of ^2'. The XIS was operated in normal 
clocking mode (8 s exposure per frame), with the stan- 
dard 5x5 and 3x3 editing mode. We processed the 
XIS data using the "xispi" and "makepi" ftool tasks and 
CALDB files of version 2009-08-13. The XIS data were 
then cleaned by assuming thresholds on the Earth eleva- 
tion angle of > 5° and the Day-Earth elevation angle of 
> 20°. We also discarded data with time since south 
Atlantic anomaly passage of less than 436 sec. After this 
screening, the remaining good exposure was 101.1 ksec 
for both FIs and BI summed over the three observations. 
Event screening with cut-off rigidity was not performed. 

The spectral analysis was performed with HEAsoft 
version 6.10 and XSPEC 12.6. To subtract the non- 
X-ray background (NXB), we employed the Dark-Earth 
database using the "xisnxbgen" ftool task (Tawa et al. 
2008). We generated two different ancillary response files 
(ARFs) for the spectrum of each region; one assumed uni- 
form sky emission, while the other utilized the observed 
XIS1 image by the "xissimarfgen" ftool task (Ishisaki et 
al. 2007). In the ARFs, we also included the effect of con- 
tamination on the optical blocking filters (OBF) of the 
XIS. 

2.2. XMM- Newton 

As mentioned in Tsuru et al. (2007), since there is a 
luminous point source within the outflow of M 82, we an- 
alyzed the XMM-Newton data to estimate its contamina- 
tion. Although XMM-Newton observed M 82 four times, 
we analyzed the data from 2004 April. This observation 
has the single largest exposure time, which is sufficient 
for this analysis. The log of this observation is given in 
table 1. The European Photon Imaging Camera (EPIC), 
aboard XMM-Newton, has two metal oxide semiconduc- 
tor (MOS) CCD arrays and one PN CCD array. All three 
cameras have a ~30' diameter circle FOV with a HPD 
of <~6", which is much smaller than the XRT and hence 
advantageous for observations of point sources. 

All data were processed using version 10.0.0 of the 
XMM-Newton Science Analysis Software (SAS). Periods 
with a 10-15 keV background count rate deviating by 
more than 3cr from the mean rate were filtered out. The 
resultant exposure times are also listed in table 1. 

3. Analysis and Results for the Disk region 

3.1. Extraction Regions of Spectra 

Figure 1 (a) shows the Suzaku X-ray image in energy 
range of 0.3-2 keV range, along with the spatial regions 
used for spectral analysis. We defined the square re- 
gion within 3'xl2' of the center of M 82 as Disk re- 
gion. Moreover, because there may be spatial variations 
in abundance patterns (Ranalli et al. 2008), the region 
between Disk and 'cap' were divided into three regions 
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Table 1. Logs of Suzaku and XMM-Newton observations 
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Instrument 


Observation 


Observation 


Effective 




ID 


Start Date 


Exposure (ks) 


XIS 


100033010 


2005/10/04 


32.3 




100033020 


2005/10/19 


40.4 




100033030 


2005/10/27 


28.4 


EPIC 


0206080101 


2004/04/21 


66.8/68.2/52.9* 



* The exposure times of the MOS1/MOS2/PN. 



as shown in figure 1 (a). These three regions are labeled 
Wind-1, Wind-2, and Wind-3 in order, from the center. 
In this paper, we derived abundance patterns of these re- 
gions. The three regions outside the M 82 wind are se- 
lected to extract background emission, and are shown as 
green ellipses (figure 1 (a)). We added the spectra of these 
three regions in order to derive the background spectra. 

In the figure 1 (b), we show the XMM-Newton EPIC 
image in 0.3-5 kcV range, overlaid with the extraction 
regions for Suzaku. As mentioned in Tsuru et al. (2007), 
the luminous X-ray source in the outflow lies in the region 
Wind-3, and we refer to it as Source A, hereafter. 

3.2. Fitting the Spectra 

We fitted the Disk region spectra with following model: 
phabsM82 x vapeciT, 2T, or 3T + phabsx-i x cutoff 
PLx-i- The phabsM82 and phabsx-i factors represent the 
column density. The first of these is left free and the sec- 
ond one is fixed to 1.1 x 10 22 cm~ 2 (Miyawaki et al. 2009). 
The ISM emission of M 82, vapecix, 2T, or 3T, was repre- 
sented by one-temperature (IT), two-temperature (2T), 
and thrcc-tcmpcrature (3T) models, employing the vapec 
plasma code (Smith et al. 2001). We divided the metals 
into the groups of O, Ne, (Mg & Al), Si, S, (Ar & Ca) and 
(Fe & Ni) based on the metal synthesis mechanism of SN, 
and allowed each group to vary. In the 2T or 3T mod- 
eling, the vapec components were assumed to share the 
same elemental abundances. Because our purpose is to 
derive metal abundance of CIE hot plasma, we accepted 
the best fit absorbed cutoff power-law (PL) model from 
Miyawaki et al. (2009) for the hard component, which is 
dominated by emission from X-l. The cutoff energy is 
fixed to 5.8 keV, which is the best-fit value in Miyawaki 
et al. (2009). Meanwhile, the photon index and normal- 
ization are set free. We simultaneously fitted the spectra 
of all CCDs, using the energy ranges of 0.4-6 keV for the 
BI CCD and 0.5-6 keV for the FI CCDs. Because there 
are known calibration issues of the Si edges in all CCDs, 
the energy range of 1.84-1.86 keV is ignored, hereafter 
(Koyama et al. 2007). 

The IT, 2T, and 3T models failed to reproduce the 
spectra, x 2 / d -°- f — 2 376/942, 1784/940, and 1673/938 re- 
spectively, despite improving fit statistic values. There 
are large residuals around 0.5-0.7 and ~1.2 keV as shown 
in top of figure 2. Although we tested fits with 4T and 
5T models, the residuals features and fit statistics do not 
change significantly, x 2 / d -°- f - =1 647/936 and 1653/934, 



respectively. Therefore, we treated thermal components 
represented by 3T models. The large residuals existing 
around 0.5-0.7 and ~1.2 keV correspond to O VII/VIII and 
NeX Ly/3 emission lines, respectively. We considered that 
these highly ionized O and Ne may occur CX interacting 
with neutrals. Adding six gaussians to represent O and Ne 
emission lines at 0.574 keV (O VII Ka), 0.66 keV (O VII K/3 
+ OVIII Lya), 0.774 kcV (OVIII Ly/3), 0.921 kcV (NelX 
Ka), ~1.02 keV (NelX Kfi + NeX Lya), and ~1.2 kcV 
(NeX Ly/3), we fitted the spectra with the above model 
while setting O and Ne abundance of vapec models to 
zero. The Gaussian line centers for only 1.02 keV and 
1.2 keV are set free, because relatively strong Fe-L line 
emission exists in this energy ranges. The phabs compo- 
nent is fixed at 29.9 x 10 20 cm~ 2 , which is value of fitting 
with 3T models. This model produced a best fit with 
X 2 /d.o.f.=1364/933 and the derived parameters are listed 
in table 2. 

Although the large residuals disappear with the above 
fit, the discrepancies of residuals between XIS0 and XIS1 
still remain, over the energy range of 0.5-0.7 keV. These 
residuals were attributed to contamination on the optical 
blocking filters (OBF) on the XIS. It is known that con- 
taminant has been accumulating on the OBF since the 
detector doors were opened following launch (Koyama et 
al. 2007). To estimate the effect of absorption, we added 
the model 'varabs', which treats photoelectric absorption 
with variable abundances. The parameters of C is set free. 
Meanwhile, we set O to 1/12 of C in units of the solar ra- 
tio, assuming the contaminant material of the OBF to be 
C24H 38 04 and owing to the solar units of the abundances 
in the model. We fitted the spectra with a model of 3T 
plus gaussians multiplied varabs, which yielded x 2 /d.o.f. 
of 1293/929. The fitted spectra are shown in the bottom 
panel of figure 2. The column density of C and O are 
2.2-7.8 and 0.37-1.33 x 10 17 cm" 2 , respectively. These 
values are consistent with those during other early obser- 
vations (Fujimoto et al. 2007). The derived parameters 
are listed in table 2. Although the fits are not formally 
acceptable, these results arc useful for assessing whether 
CX is important or not. 

4. Analysis and Results for the Wind regions 

4-1. Contamination of Disk Emission 

The Disk emission contaminates the Wind regions ow- 
ing to the low Suzaku angular resolution (HPD ~ 2'). We 
simulated the emission from the Disk region in each Wind 
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Fig. 1. Panel (a) and (b) are the X-ray images of M 82 with the Suzaku XIS1 and XMM-Newton EPIC-MOS1, respec- 
tively. The background components were not subtracted and vignetting was not corrected for in these images. Panel (a); 
the XIS 1 image in 0.3-2 keV is smoothed with a <r=5.2" Gaussian profile and overlaid with optical galactic disk con- 
tours from the Digitized Sky Survey (yellow). The magenta boxes show the regions for spectra extraction labeled Disk, 
Wind-1, Wind-2, and Wind-3 going out from the center. The background spectra were extracted from the three green 
ellipses. The gray ellipse shows the 'cap' region (Tsuru et al. 2007). Panel (b); the selected Wind regions are over- 
laid. The position of Source A is also indicated. The white box is the Suzaku XIS field of view in both panels. 



region, setting the point source to the coordinates of X-l 
by xissimarfgen in the energy range of 0.2-16 keV, because 
the Suzaku PSF does not allow us to resolve the X-ray 
peak of hot ISM and X-l (Ishisaki et al. 2007). We found 
a level of 9%, 1%, and 0.5% Disk intensity contamination 
to the Wind-1, Wind-2, and Wind-3 regions, respectively. 
The level of contamination is almost constant, to within 
5% below 6 keV, and within 20% at 6-10 keV. As shown in 
figure 3, most of the Wind-1 emission beyond 1.5 keV can 
be attributed to the Disk region. Meanwhile, in the Wind- 
2 and Wind-3 regions, the contamination from Disk emis- 
sion is smaller compared with Wind-1. The intensity of 
contamination from Disk region is one-half and one- fifth of 
that in Wind-2 and Wind-3 above ^1.5 keV, respectively. 
Hereafter, we analyze the spectra of Wind-1 and Wind- 
2 regions after subtraction of the directly-scaled emission 
of the Disk region. Full statistical error propagation was 
employed in this analysis. 

4-2. Wind-1 Region 

We first fitted the Wind-1 region spectra after subtract- 
ing 9% of the Disk emission and the background emission 
accumulated over three background regions with the fol- 
lowing model; phabsc X vapec2T- The Wind-1 region 
spectra will be reproduced as the sum of Disk emission, 
power-law component of X-l, and Wind-1 emission. The 
phabsc factor is fixed to 4.0 x 10 20 cm~ 2 , which is the 
Galactic value of TVh in the direction of M 82. The abun- 
dances of O, Ne, Mg, and Fe are allowed to vary, and 
those of other metals are fixed to solar value. We simul- 



taneously fitted the spectra of all CCDs over the energy 
ranges of 0.4-5.0 keV for the BI CCD and 0.5-5.0 keV 
for the FI CCDs. The IT and 2T models failed to repro- 
duce the spectra, with x 2 /d-o.f.=845/305 and 807/303, 
respectively. Above 2 keV, there are large negative resid- 
uals which we attribute to excess subtraction of the Disk 
emission contamination. 

The simulated intensity of Disk contamination may 
have some uncertainties due to spatial spread of the X-ray 
emission and uncertainties in the point spread function. 
Therefore, we tested the effect of varying the Disk contam- 
ination level over 6-8.5%, subtracting off this component, 
and then carrying out the fit. The x 2 /d.o.f. values are 
plotted in figure 4, and are found to show a minimum 
at 354/294 with the contaminant level of 7% of the Disk 
emission. Above 2 keV, the spectra of the Wind-1 region 
are well reproduced by the contamination of 7% emis- 
sion from the Disk region (figure 3), and these parameters 
were adopted hereafter. The derived parameters and fit- 
ted spectra are shown in the first row of table 3 and in 
figure 5, respectively. As shown in figure 5, in the Wind- 
1 spectra, several emission lines are detected. The ones 
around 0.5-0.6 keV, 0.6-0.7 keV and ~1.3 keV are identi- 
fied with Ka lines of O VII, O VIII, and Mg XI, respectively. 
The emission bump around 0.7-1 keV corresponds to Fe-L 
complex, as well as to K-lines from Ne IX and Ne X. 

To investigate the sensitivity of metal abundance deter- 
mination on the ISM temperature, we also employed the 
IT and 3T models and fitted the spectrum again, assum- 
ing that the Disk contamination level is 7%. The IT and 
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Fig. 2. Background-subtracted XISO (black) and XIS1 (red) spectra of Disk, shown without removal of the instrumen- 
tal response. In order from top to bottom, the panels are the spectra fitted with model 3T for ISM, 3T + Gaussian 
lines, and 3T + Gaussian lines with absorption of OBF contaminant. Black and red lines show the best-fit model 
in each case for FIs and XIS1, respectively. For simplicity, only the model components plot of the XIS1 spectra arc 
shown. Light blue, blue and purple lines are the three ISM temperature components. Green curves represent the accu- 
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Fig. 3. From top to bottom, panels are the comparison spectra between scaled contamination Disk and Wind-1, 
Wind-2, and Wind-3, respectively. In the top panel, the spectra of Wind-1 and Disk scaled by 9% and 
7% are shown in black, light blue, and blue, respectively. The contamination-subtracted Wind-1 and back- 
ground region spectra are shown as red and green, respectively. The other panels follow the same layout. 
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Table 2. Summary of the best-fit parameters for the Disk region with 3T + Gaussian models. 
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P a r m o t, o r s 




Disk 
3T for ISM 


Gaussians 


Gaussians 
(contami) 


Nvt 

± * il]VI82 


(10 20 cm- 2 ) 


— 1 20 


29.9 (fix) 


29.9 (fix) 




(keV) 


1 O3+0:n 


1 7 n+0.03 


1 fi2+ 004 

— 08 


kT 2T 


(keV) 


0.59±0.01 


0.61±0'.01 


0.58±0.01 


kT 3T 


(keV) 


0.24±0.02 


n 34+00 4 

u -°^-0.02 


0.08±0.01 





(solar) 


n 4^+0-05 

u - 4o -0.04 


(fix) 


(fix) 


Ne 


(solar) 


1 iq+0. 09 
i - 1J -0.08 


(fix) 


(fix) 


Mg, Al 


(solar) 


1.04±0.08 


1 10+0. 09 
i - lo -0.05 


1 - Ua -0.11 


Si 


(solar) 


1 9 o+o.o9 

1 - ZJ -0.08 


1 9 o+0.09 


1 10+0. 03 
i - i8 -0.05 


S 


(solar) 


1.52±0.11 


1 4n +0 - 08 

1 - 4U -0.08 


1 OO+0.07 
1 -°°-0.08 


Ar, Ca 


(solar) 


0.65±0.31 


1 -ifi+o- 27 

l.OD_ 26 


1 OO+0.32 
l.OO_ 26 


Fe, Ni 


(solar) 


0.38±0.03 


n 42+ - 05 

U-^-o.Ol 


n 42+ - 03 



Tcutoff PL 

line energy 
Norm 
line energy 
Norm 
line energy 
Norm 
line energy 
Norm 
line energy 
Norm 
line energy 
Norm 



0.47 



+0.01 

-0.03 



0.52±0.03 0.55 



+0.02 



(keV) 

(10- 4 ) 
(keV) 

(10- 4 ) 
(keV) 

(10- 4 ) 
(keV) 

(10- 4 ) 
(keV) 

(10- 4 ) 
(keV) 

(10- 4 ) 



0.574 (fix) 
I.28±0.16 
0.66 (fix) 
3.34±0.13 
0.774 (fix) 
I.08±0.12 
0.921 (fix) 
2.24±0.10 
1.03±0.01 
3.20±0.09 
1.24±0.01 
0.95±0.05 



0.04 

0.574 (fix) 

1 or+0.16 

1 - oo -0.15 

0.66 (fix) 

o 4 r;+0.14 

0.774 (fix) 

1.16±0.12 

0.921 (fix) 

6 <S2 +0 - 20 
°- oz -o.io 

1.02±0.01 
3.22±0.09 
1.24±0.01 
0.99±0.05 



X 2 /d.o.f. 



1673/938 1364/933 1293/929 



3T model fits result in x 2 /d.o.f. =479/296 and 351/292, 
respectively. In particular, IT model gives a much worse 
fit to the data as compared to the 2T model which had a 
X 2 /d.o.f.=354/294. The derived parameters are summa- 
rized in table 3. Furthermore, we calculated confidence 
contours between the metal abundances (O, Ne, and Mg) 
relative to Fe, using the IT, 2T, and 3T ISM models. 
From these contours, we derived the abundance patterns 
(table 3), which have smaller uncertainties as compared 
to the absolute abundance values. Figure 6 shows the 
abundance patterns of M 82 Wind-1 with IT, 2T, and 3T 
models. These metal abundance patterns are consistent 
within error in all models, although the 3T model shows 
the largest uncertainties. 

4.3. Wind- 2 Region 

We fitted the Wind-2 region spectra after subtracting 
1% of the Disk emission and the background emission ac- 
cumulated over three background regions with the same 
model as the Wind-1 region. In this case, the energy 
ranges are 0.4-5.0 keV for the BI CCD and 0.5-5.0 keV 
for the FI CCDs. The fit statistic of the 2T model was 
not good at x 2 /d.o.f.=397/292. Because positive residuals 
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Fig. 4. The x 2 /d.o.f. 
when fitting the Wind-1 
traction of various Disk 



distribution obtained 
spectra, after sub- 
contamination levels. 
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Fig. 5. Background-subtracted XISO (black) and XIS1 (red) spectra of Wind-1 (left) and Wind-2 (right), shown without re- 
moval of instrumental response. Black and red lines show the best-fit model for the XISO and XIS1, respectively. For 
simplicity, only the model components for XIS1 spectra are shown. Blue and light blue lines are the ISM components. 



Table 3. Summary of the best-fit parameters for the Wind-1 region with IT, 2T, and 3T for ISM models. 



Parameters 



Wind-1 
IT for ISM 
Disk 7% 



2T for ISM 3T for ISM 



(keV) 
(keV) 
(keV) 
(solar) 
(solar) 
(solar) 
(solar) 
(solar) 
(solar) 
(solar) 



kT 1T 

kT 2 T 

kT 3T 

O 

Ne 

Mg, Al 
Fe, Ni 
O/Fe 
Ne/Fe 
Mg/Fe 

X 2 /d.o.f. 



4.0 (fix) 
0.38±0.01 



5+0.10 
-0.09 
1+0.14 
-0.13 

o.85±S:iS 

0.35±0.04 

9 r+0.33 
z -°-0.25 
+0.21 
-0.36 
-0.28 
-0.42 



0. 
1.29! 



3.86 
2.5" 



4.0 (fix) 
0.26±0.01 
0.57±0.02 



+0.17 
-0.12 
-0.29 
-0.22 
0.23 
-0.18 
0.08 



1.00 
1.73 
1.09 
0.51 
1.96 ' 
3.39 
2.13 



0.25 

+0.54 
-0.50 
-0.37 
-0.30 



4.0 (fix) 
0.26±0.01 
0.54^ 
0.97^ 
0.99^ 
1.75" 1 



,+0.04 
-0.05 
-1.04 
-0.21 
1-0.16 
-0.13 
h0.32 

'-0.26 
i 9 r+0.32 
1 - ZO -0.24 
1-0.14 
-0.12 
+ 1.01 
0.40 
+0.74 
-0.43 
-0.34 
-0.33 



0.60 
1.65 
2.93 
2.08 



479/296 354/294 



351/292 



exist above 2 keV, we attempted setting the contamina- 
tion at 1.5% and 2% intensity of the Disk emission. After 
subtracting and refitting, we found x 2 /d.o.i. of 360/303 
and 668/295, respectively. We accepted 1.5% as the con- 
taminant ratio of Disk emission to Wind-2 region. The 
resultant parameters and fitted spectra are summarized 
in table 4 and figure 5, respectively. 

In Wind-2 region, we also employed the IT and 3T 
models in fitting, assuming that the Disk contamination 
level is 1.5%. The fitting of IT and 3T model gives 
X 2 / d -°- f -= 448 / 305 and 352/301, respectively. Again, the 
IT model returns a worse x 2 /d.o.i. of the 2T case 
360/303. The metal abundance ratios are derived from 
confidence contours and are plotted in figure 6. Although 
the O/Fe ratio of the 3T model is slightly lower than the 
others, there is no large discrepancy amongst the various 



temperature models. 
4-4- Wind-3 Region 

4- 4-1. XMM-Newton Analysis of Source A Spectra 

Before fitting the Wind-3 spectra, we estimate the emis- 
sion from Source A based upon XMM-Newton data. We 
extracted source spectra within a 15" radius circular aper- 
ture centered at (a, 5) = (09 h 55 m 14 s , +69°47'35"). The 
background spectra were extracted from an annular ring 
and subtracted from Source A. We extracted the spectra 
of MOS 1, MOS 2, and PN and fitted these simultane- 
ously in the 0.3-7.0 keV energy range with the following 
model; phabsA x PLa, with all parameters set free. The 
best fit values are listed in table 5. The Source A spectra 
are reproduced by the absorbed PL model with a photon 
index of l-76±g;Jg, which is consistent with a typical active 
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Fig. 6. Abundance ratios of O, Ne, and Mg to Fe for the IT, 2T, and 3T model of the M 82 Wind-l 
(left) and Wind-2 (right) regions. The abundance patterns of Wind-l and Wind-2 region are also shown. 
Solid, dot-dashed, and dotted lines represent the number ratios of metals to Fe for the solar abundance, 
for SN II, and for SN la products (Lodders 2003; Iwamoto et al. 1999; Nomoto et al. 2006), respectively. 



Table 4. Summary of the best-fit parameters for the Wind-2 region with IT, 2T, and 3T for ISM models. 



Parameters Wind-2 

IT for ISM 2T for ISM 3T for ISM 
Disk 1.5% 





(10^ u cm" 2 ) 


4.0 (fix) 


4.0 (fix) 


4.0 (fix) 


kT 1T 


(keV) 


0.32±0.01 


0.24±0.01 


0.24±0.01 


kT 2 r 


(keV) 




n co+0.06 

u.oa_ 05 


™+ - 08 


kT 3 T 


(keV) 






1 1S+0- 47 


O 


(solar) 


61+ 09 

u - oi -0.07 


69+ 012 


69+ 013 


Ne 


(solar) 


u ' y ^-0.12 


1 20+ 23 

x - zu -0.19 


-, oc+0.31 
L - Zb -0.25 


Mg, Al 


(solar) 


70+ - 20 

u - ' u -0.17 


62+ - 21 

u - oz -0.17 


0.84lg;p 


Fe, Ni 


(solar) 


u -°^-0.05 


n OO+0.08 
U-JO-O.Og 


^7+°- 21 

u -°'-0.15 


O/Fe 


(solar) 


1 7fi+0- 24 

x - '°-0.26 


i o 9 +0.40 


9 +0.51 


Ne/Fe 


(solar) 


z -°' -0.23 


q 1 O + 0.57 

°- lo -0.84 


9 99+0.64 
z - zz -0.72 


Mg/Fe 


(solar) 


i 4 +0.36 


1 61 + - 47 

1 - D1 -0.50 


1 4t;+0.55 



X 2 /d.o.f. 



448/305 



360/303 



352/301 
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Table 5. Best-fit parameters for Source A with phabsA X 
PL A 



Parameters 


^H A 

r A 


(10 20 cm" 2 ) 


i i 9 +2. 2U 
11>z -2.44 

1 7fi+ 010 

i -'°-0.09 


Source A flux* 


(erg cm -2 s -1 ) 


2.2t^xl0- 13 


X 2 /d.o.f. 




167/178 



* Flux within the accumulated region between 0.5 and 8 keV. 



galactic nucleus spectrum. 

4-4-2. Estimation of Background Emission 

As a result of lower surface brightness of the Wind- 
3 region as compared to Wind-1 and Wind-2 (figure 3), 
we decided to perform simultaneous fitting of the Wind-3 
and background spectra. Background spectra were ex- 
tracted from three regions. Because the NXB is sub- 
tracted, the background components consist of Cosmic 
X-ray Background (CXB) and Galactic thermal emission. 
We represented CXB with a PL model and Galactic ther- 
mal emission as a two-temperature apec model (Smith et 
al. 2001). Empirically, one thermal component represents 
the sum of solar wind charge exchange (SWCX) and local 
hot bubble (LHB), while Milky Way halo (MWH) emits 
the other thermal component (Yoshino et al. 2009). 

We fitted the spectra with the following models; phabsc 
x (PLcxb + apecMWH) + apecLHB • The 'phabsc' model 
represents the photoelectric absorption, whose column 
density was fixed to the Galactic value of 4.0 x 10 20 ct7j- 2 
in the direction of M 82. The 'apecMWH' and 'apecLHB' 
models correspond to thermal emission from our Galaxy, 
with the metal abundances fixed at the solar value and 
zero redshift. The spectra from the BI and FI CCDs 
were simultaneously fitted in the 0.4-5.0 and 0.5-5.0 keV 
range, respectively. In table 6, we list the derived pa- 
rameters, which are consistent with Tsuru et al. (2007). 
This model reproduced the spectra approximately, with 
X 2 /d.o.f.=312/255. The derived photon index of PL CXB 
is 1.4310°^, which is consistent with Kushino et al. (2002). 
Although the temperature of apecMWH model is higher 
than that of the common MWH component (^0.25 keV), 
there are known to be high temperature Galactic com- 
ponent (~0.6 keV) present in some portions of the sky 
(Yoshino et al. 2009). Furthermore, similar results are 
reported in Sato et al. (2008) and Konami et al. (2010). 
Since the Galactic emission is considered to have spatial 
dependencies, we regard this model to be appropriate for 
estimating the background emission in this analysis. 
4-4-3- Simultaneously Fitting of Wind-3 and Background 
Regions 

We fitted the Wind-3 spectra and background simulta- 
neously with following models; phabsc x vapecix,2T + 
phabsA x PLa + background models. Here, we used en- 
ergies ranging from 0.4 (BI) or 0.5 (FI) to 5.0 and 8.0 keV 
for the background and Wind-3 regions, respectively. The 
first and second terms represent the outflow and Source- 
A, respectively. The photon index of the PL model for 



[Vol. , 

Table 6. Best-fit parameters for the background regions 
with apec components + PL model.* 



Parameters 


Tcxb 


(10^ u cm" 2 ) 


4.0 (fix) 

1 43+009 


kTyiWH 

Abundance 

fcTLHB 

Abundance 
A r orm LH B/A f orm M wH 


(keV) 
(solar) 
(keV) 
(solar) 


0.70±0.07 

1 (fix) 
21 + 03 

u - zl -0.02 

1 (fix) 

0.97i°;i 2 


X 2 /d.o.f. 




312/255 



* The apec components for spectra in the background region of 
M 82 with absorbed MWH component, LHB component for 
the Galactic emission, and a PL model for CXB. 



CXB was fixed to 1.43, which was best fit parameter in 
section 4.4.2. All the other background parameters were 
linked between the two regions. 

We employed the IT and 2T models in fitting, and 
found x 2 /d.o.f. = 601/543 and 587/541, respectively. The 
fitted spectra and derived parameters are shown in figure 
7 and table 7, respectively. The temperatures and metal 
abundances between the IT and 2T models are consistent 
with each other, although the 2T model allows higher val- 
ues of abundances. When employing the 3T model, the 
third temperature could not be constrained and the nor- 
malization is less than 10% of the those of other compo- 
nents. The spectra of Wind-3 region do not require a third 
temperature component. 

4-5. Search for CX emission in Wind regions 

As seen in section 3.2, the CX emission may contribute 
in addition to thermal emission in Disk region. To search 
for CX emission in Wind regions, we fitted the spectra 
of Wind-1, Wind-2, and Wind-3 regions with a 2T ther- 
mal model including six Gaussians representing O and Ne 
emission lines, by setting O and Ne abundance of vapec 
models to zero, as described in section 3.2. The ratios of 
(OVII K[3 + OVIII Lya)/0 VII Ka, OVIII Ly/3/(OviI K>3 
+ OVIII Lya), (NelX K/3 + NeX Lya)/NeIX Ka, and 
NeX Ly/3/(NeIX K/3 + NeX Lya) are derived in Disk, 
Wind-1, Wind-2, and Wind-3 regions. 

5. Discussion 

5.1. Excess emission of Lyj3 lines of H-like O and Ne 

There are residual structures at 0.5-0.7 keV and ~ 1.2 
keV in the spectra of the Disk region fitted with the ther- 
mal models. These residuals are well fitted by excess emis- 
sion of Ly/3 lines of O VIII and Ne X. A plausible expla- 
nation for these residuals is CX between highly ionized 
O/Ne and neutral H/He, since the CX process gives higher 
ratios of K/3/Ka or K7/K/3 than those of CIE plasma. 
We calculated confidence contours between the normal- 
izations of the Gaussian lines, OVIII Ly/3/(OviI K/3 + 
OVIII Lya) and NeX Ly/3/(NeLX K/3 + NeX Lya), us- 
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Fig. 7. NXB-subtracted XISO (black) and XIS1 (red) spectra of M82 (left panel), and those of the background region (right panel), 
shown without removal of instrumental response. The left panel employs the two-temperature model for the ISM. Black and red lines 
show the best-fit model for XISO and XIS1, respectively. For simplicity, only the model components for XIS1 spectra are shown. Blue 
and light blue lines are the ISM components, the green line is emission from Source A modelled as an absorbed power-law (PLa), 
gray lines are the Galactic background emission (apec^wH an d apecLHB) an d the CXB components, respectively. The background 
components are common between the on-source and background spectra, but scaled to the respective data accumulation area. 



Table 7. Summary of the best-fit parameters for the Wind-3 region with IT or 2T for ISM models. 



Parameters 




Wind-3 




AT 


Cir.aU „™-2n 


IT for ISM 


2T for ISM 



He 

kT 2 r 

Normn /Norrri2T 

O 

Ne 

Mg, Al 
Fe, Ni 
O/Fe 
Ne/Fe 
Mg/Fe 

^H A 

r A 

Source A flux* 



cm 
(keV) 
(keV) 

(solar) 
(solar) 
(solar) 
(solar) 
(solar) 
(solar) 
(solar) 

(It) 20 cm^ 2 ) 
(erg cm~ 2 s - 



n 41+0-05 

-0.03 



0.57± ; 2 ? 
61 +0 - 25 

u - zl -0.07 
, + 1.66 

-0.98 



2.73] 
4.55 
2.90 



) 2.1 



11.2 
1.73 

h0.5 
-0.2 



+ 1.70 

-1.30 
+ 1.39 
-1.15 

(fix) 

+0.23 
-0.16 



20+ - 04 

u ' zu -0.06 

4Q+0-08 

u -^-0.09 

U '°°-0.35 
+0.55 
-0.26 
1.21 
0.47 
+0.64 
-0.34 

u - ol -0.07 
7+0.77 

-0.97 

+ 1.10 
-1.15 



0.68 
1.48 
0.90 



2.17] 



2.90 



xlO 



-13 



11.2 (fix) 

1.6818:3 
2.i^: 6 39 xio- 13 



X 2 /d.o.f. 



601/543 



587/541 



* Flux within the accumulated region between 0.5 and 8 keV. 
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ing thermal models including six Gaussians representing 
O and Ne emission lines (see subsection 3.2). In figure 8, 
the derived these values are plotted against plasma tem- 
perature with those of the CIE plasma from the APEC 
code. Both these ratios for H-like O and Ne in the Disk 
region have values of ~0.3, significantly larger than those 
predicted for the CIE plasma. Since the temperature de- 
pendence of Ly/3 to Lya ratio is relatively small, no com- 
bination of temperature components alone can give the 
observed high Ly/3 excess. In contrast, the Wind regions 
have ratios of OVIII Ly/3/ (O VII K/3 + OVIII Lya) and 
NeX Ly/3/ (Ne IX K/3 + NeX Lya) which are consistent 
with those for the CIE plasma. Especially, in the Wind-1 
region, OVIII Ly/3/(OviI K/3 + OVIII Lya) is factor of 
3-10 smaller than those in the Disk region. Therefore, 
this line ratio analysis indicates that the CX process gives 
higher Ly/3 to Lya ratio of H-like O and Ne in the Disk 
region. 

For He-like ions, the values of (Ovil K/3 + OVIII 
Lya)/0VII Ka and (NelX K/3 + NeX Lya)/NeIX Ka 
in the Disk and the all Wind regions are consistent with 
those of CIE plasma. The cross-sections of interaction of 
highly ionized H-like ions with neutral H/He atoms are 
larger than those of He-like ions (Koutroumpa et al. 2006 
and references therein). Since K/3 line energies of He- like 
O and Ne are close to those of Ka lines of H-like ions of the 
same element, and since the ratios of He-like to H-like ion 
Ka lines have very steep temperature dependencies, K/3 
to Ka ratios of He-like ions are not suitable for the study 
of the contribution of CX emission using CCD spectra. 

The Ly/3 lines of H-like O and Ne lie close the expected 
energies of strong Fe-L lines. Large residuals at ~ 1.2 
keV have been observed in various kinds of objects (e.g. 
Yamaguchi et al. 2010). In elliptical galaxies with hot 
ISM of ~ 0.6 keV, residuals at ~ 0.8 keV are also observed 
(Matsushita et al. 2007; Tawaraet al. 2008). These residu- 
als have been interpreted as a problem in the Fe-L atomic 
data (e.g. Brickhousc et al. 2000). However, there is a 
clear difference in the residual structures between the Disk 
and the Wind-1 regions of M 82. In these two regions, the 
O and Ne lines mostly come from the same two temper- 
ature components, 0.2-0.3 keV and ~ 0.6 keV, although 
the ratio of the normalization is different. Therefore, this 
spatial variation of the residual structure indicates the dif- 
ference to be a result of CX rather than uncertainties in 
the Fe-L atomic data. 

Liu et al. (2011) discovered excess emission of the for- 
bidden lines in the triplet of the He-like O and Ne, which 
are robust evidence of CX. They assumed the temperature 
of thermal emission is ~0.6 keV, which is derived using the 
ratio of MgXII Lya/MgXI Ka. Therefore, the CX emis- 
sion of Ovil is a result of the capture of two electrons 
by O IX ions, while the cross-section of this process is less 
than 20% of that of single-electron capture via the CX pro- 
cess of O VIII, because the ion fraction of O IX is dominant 
at the ~0.6 keV CIE plasma (Liu et al. 2011, Greenwood 
et al. 2001). Liu et al. (2011) evaluated the contribution of 
CX to be ~90% to the O VII triplet, assuming the temper- 
ature of thermal emission to be ^0.4-0.6 keV. The contri- 



bution of CX depends on the temperature structure, since 
the fraction of resonance, intercombination, and forbidden 
lines are dependent on the plasma temperature (Porquet 
& Dubau 2000). The ISM of galaxies is empirically as- 
sumed to be composed of a multiple temperature plasma. 
In fact, our results need ~0.2 keV temperature plasma for 
the disk emission, which accounts for ^30% of the He-like 
triplet of O. The existence of the 0.2 keV plasma increases 
the ratio of the forbidden to resonance lines, but this is 
still much smaller than the observed values in Liu et al. 
(2011). With Astro-H, next Japanese X-ray satellite, we 
will be able to resolve these Ka and K/3 lines of H-like 
and He-like ions, as well as the triplets of He-like lines, 
and investigate CX more clearly. 

CX should strongly contribute to the X-ray emission in 
the central regions of starburst galaxies, where abundant 
cold gas and highly ionized ions co-exist and interact. We 
have detected plausible evidence of CX emission between 
H-like ions (O and Ne) and neutrals (H and He) in the 
Disk region. In the 'cap' region, the ionized outflow from 
M 82 are thought to collide with cool ambient gas and 
enhance the Ha and X-ray emission, and may also cause 
CX emission (e.g. Devine & Bally 1999, Lehnert 1999). 
An emission line consistent with the C VI transition of n = 
4 to 1 at 0.459 keV has been marginally detected, although 
it is not statistically significant at the 99% confidence level 
(Tsuru et al. 2007). 

5.2. Metal Abundance Patterns of All Wind regions 

We have successfully derived metal abundances of O, 
Ne, Mg, and Fe over the entire Wind region of M 82 for 
the first time. The spectra are well reproduced with ther- 
mal models and without evidence of CX emission in the 
outflow regions. In these regions, the amount of neutral 
gas may be smaller than those in the Disk and 'cap' re- 
gions, and therefore, the effect of the CX may not be im- 
portant. As a result, the metal abundances in the Wind 
regions should be more reliable than those in the Disk and 
'cap' regions. Furthermore, there no discrepancies among 
the metal abundance ratios determined using either the 
IT, 2T, or 3T ISM model over all the Wind regions. We 
calculated confidence contours between the abundance of 
metals (O, Ne, and Mg) relative to Fe, using the 2T model 
for ISM of all Wind regions. From these contours, we de- 
rived the abundance patterns (table 3, table 4, and table 
7), which have smaller uncertainties as compared to the 
absolute abundance values. Figure 9 shows the abundance 
patterns of M 82 Wind-1, Wind-2, and Wind-3 within the 
'cap' region (Tsuru et al. 2007). The distance of these re- 
gions from the galactic center range from 2 kpc to 10 kpc. 
The O/Fe, Ne/Fe, and Mg/Fe ratios are consistent with 
no spatial variation. 

The calculated SN II and SN la yields are also plotted in 
figure 9. The SN II yields by Nomoto et al. (2006) refer to 
an average over the Salpeter initial mass function of stellar 
mass from 10 to 50 M©, with a progenitor metallicity of 
Z = 0.02. The SN la yields were taken from the W7 model 
Iwamoto et al. (1999). The O/Fe, Ne/Fe, and Mg/Fe ra- 
tios are a factor of 2-3 higher than the solar values, and 
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Fig. 8. Left; The red, green, blue, and light blue data represent the values of the line intensity ratio Oviii Ly/3/(Ovn K/3 + 
O VIII Lye?) vs. the thermal model temperatures, for the Disk, Wind-1, Wind-2, and Wind-3 regions, respectively. The temperatures 
used are those determined from thermal models including six Gaussians representing O and Ne emission lines, by setting O and 
Ne abundance of vapec models to zero. The solid line is same ratio predicted for a CIE plasma using the vapec code (Smith 
et al. 2001). Right; Symbols and curves are as in the left-hand panel, but for the line ratio Nex Ly/3/(Nelx K/3 + Nex Lya). 



therefore, the outflow of M 82 is significant enriched by SN 
II yields. The outflow wind of another starburst galaxy, 
NGC 4631, observed with Suzaku also shows enhance- 
ments of a-elements relative to the Fe abundance ratios 
(Yamasaki et al. 2009), although its disk region possesses a 
hot ISM with an approximately solar abundance pattern. 
The abundance pattern of the hot ISM in NGC 4258, a 
spiral galaxy without starburst activity, is also consistent 
with the solar ratio (Konami et al. 2009). These results 
indicate that starburst galaxies eject metals synthesized 
by SN II into the intergalactic medium via outflows. The 
derived abundance pattern of the Disk region, especially 
the very low O abundance relative to the other elements, 
is similar to those observed with ASCA and XMM (Tsuru 
et al. 1997; Ranalli et al. 2008). However, the observed 
excess emission of Ly/3 lines and forbidden lines in He- like 
triplet indicates that a significant fraction of these lines 
are caused by CX process. Therefore, it is very difficult to 
isolate the thermal emission owing to the multiple temper- 
ature ISM plasma. Furthermore, the derived value of Nu, 
~ 30 x 10 20 cm~ 2 , is significantly higher than the Galactic 
value of ~ 4 x 10 20 cm~ 2 , whereas those of the Wind re- 
gions are consistent to the Galactic value. The observed 
value of Nh toward the Disk region means that 85% of 
photons in the Lya line of H-like O must be absorbed, 
while only 20-30% of K lines of Mg absorbed. The cold gas 
which absorbs the X-ray emission from the Disk should be 
non-uniform, and X-ray photons should be partially ab- 
sorbed. Therefore, the observed very low abundance of O 
may be caused by incorrect modelling of spectral compo- 
nents of this very complicated plasma in the Disk region. 



6. Conclusion 

We have performed X-ray spectral analysis of the out- 
flow region of M 82 observed with Suzaku. Suzaku XIS 
has a high spectral sensitivity to faint O lines thanks to a 
very small low-pulse-height tail below ~1 keV. The spec- 
tra from the galactic disk are reproduced with 3T ther- 
mal models and excess emission of Ly/3 lines of O VIII and 
Ne X, which may be caused by CX process. In the outflow 
region, there is no excess Ly/3 emission, and the derived 
O/Fe, Ne/Fe, and Mg/Fe ratios are a factor of 2-3 higher 
than the solar ratio. This abundance pattern indicates 
that starburst activity enriches the intergalactic medium 
with SN II yields via outflows. 

We thank the referee for providing valuable comments. 
We gratefully acknowledge all members of the Suzaku 
hardware and software teams and the Science Working 
Group. SK is supported by JSPS Research Fellowship for 
Young Scientists. 
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Fig. 9. Abundance ratios of O, Ne, and Mg to Fe for the best-fit model of the M 82 Wind regions. The abun- 
dance patterns of Wind-1 (green), Wind-2 (blue), Wind-3 (light blue), and cap (magenta; Tsuru et al. 2007) regions 
are also shown. Solid, dot-dashed, and dotted lines represent the number ratios of metals to Fe for the solar abun- 
dance, for SN II, and for SN la products (Lodders 2003; Iwamoto et al. 1999; Nomoto et al. 2006), respectively. 
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